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2 

A “ Cross sectional area for flow, ft 

2 

Aj - Heat transfer area, ft 

Op Specific heat at constant pressure, Btu/ltra 

Dg - Equivalent diameter, ft 

d ~ Spacing of the tubes, inch, 

dg ~ External diameter of the tubes of primary circuit, 
inch. 

d^ - Internal diameter of the tubes of primary circuit, 
inch . 

e - Energy per unit mass, Btu/lbn 

Eg - Frictional force, Ibf 

f - Friction coefficient 

G- - Mass flux, Ibm/ft^-hr 

2 

g - Acceleration due to gravity, ft/sec 

2 

g - Conversion ratio, Ibra-f t/lbf~sec 

H - Specific enhalpy, Btu/lbm 

^fg ~ batent heat of vaporization, Btu/lbm 

H^n “* Enthalpy at the inlet of the generator, Btu/lbm 

h ~ Mechanical equivalent of heat, J » 778 Ibf-ft/Btu 

k - Thermal Conductivity, Btu/hr-f 

- Length of the evaporator, ft 



P 

AP 

Q 


T 

"p 
A T 

u 

U 

V 


sat 


V 

V 

x,X 

^tt 

Z 

A Z 
a 

f 

a 


Length, of the preheater, ft 

Wetted perimeter, ft 

Pressure, Psi 

Pressure drop, Psi 

Rate of heat transfer, Btu/hr 

Heat flux, Btu/hr-ft^ 

T emp er atur e ^ 

Temperature of the wall of the tube 


- ^Sat' 

Internal energy, Btu/lhm 
Total internal energy, Btu/lbni 
Velocity, ft/hr 
Specific voliAme, ft^/lbm 
Mass flow rate, Ibm/hr 
Quality of vapour 

Locldiart and Martinelli's pareaneter 


Direction of flovf 

Height of each control volume, ft. 

Void fraction 
Viscosity, Ibf/iir-ft 
Density, Ibra/ft^ 

Surface tension, Ibf/ft 
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^£o - Multiplicator for calculating frictional 

pressure drop in two-phase flow 

Y - Empirical constant used in Thom’s correlation 

for calculating the void-fraction 


Sun erscrj-pt ; 


Refers to the value at the middle of a control 
volume. 


Subsorit>t s 

e - 

ev 

f 

g 

i 

iso 

3 

1 

p 

PA 

S 

Sat - 

V - 

c.v. 


Refers to the exit 
Evaporator 

Saturated liquid condition 
Saturated vapour condition 
Inlet condition 

Refers to isothermal conditioii 

Refers to junction 

Saturated liquid 
Refers to primary circuit 

Refers to preheater 

Refers to secondary circuit 

Saturation condition 

Saturated vapour condition 
Control volume. 



0,1 


Outlet wall condition of tlie control volume, i 


D iro en s i cn less Gr ou -p o ; 

Nu ♦” Nusselt Murabor = 1 

Pr “■ Prandtl Number = 

Re “ Reynolds Number = 


i.D^/k 
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ABSTHACT 

A digital computer code was developed by Kar [l] 
to simulate the transient behaviour of a steam generator 
typica3.1y used in Pressurized V/ater reactor type of Nuclear 
Power Plants. The present work is an extension to the above 
thesis. 

The three fundamental time dependent equations of 
conservation of mass, energy and momentum were converted into 
a set of linear algebraic equations using finite difference 
method, from which an explicit numerical scheme was formulated. 

On the primary side single phase flow and forced 
convection heat transfer were assumed and on the secondary 
side the possibility of the existance of different flow 
regimes was assumed and the heat transfer in these regimes 
was calculated using appropriate correlations. 

Some important modifications were made and the 
numerical instabilities in the program were corrected, 
program was studied for disturbances in flow and internal 
energy of varying step increments. Suitable time step was 
found. Some ranges for the step increments were fixed. 
Comparison of the results was made v^ith the two previous 
works and the results are in agreement. 
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CHAP TEH I 


1.1 INTHODUGTIOH 

V/ith tile increasing trend in the studies of 
safety measures of nuclear reactor 'beha‘viour,it becomos nneessary 
to simulate the behaviow of various components of tiie 
reactor. Dynamic analysis is an essential part of system 
design for nuclear power plants. SimuJ-ation studies prcvide 
information needed to assess plant maneuvering capability, 
control strategics and tolerance to malfunctions. 

One of the iraportant components of the nuclear 
reactor system is the steam generator, which serves as heat 
transfer link between the coolant, which flovrs througli the 
primary side and the secondary fluid, which provides power 
for the turbine. 

It is quite difficult to model a large system via 
state variable methods, because of the fact, that the overall 
system has to be broken into various sub-'Oy stems, and then 
the individual systems hai^e to be modelled seperately. The 
overall system is then constructed by coupling the sub-syston 
models. The steam generator presents to a complex modelling 
problem, as it involves sub-cooled aud boiling heat transfer 
on the secondary side. 
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DynEimic modelling is particularly important vrhen 
system design with new features are contemplated. It is 
not possible to assess the performance of the new system 
design with the previous experience. 

Though a great deal of care. is taken to ensure 
public safety in nuclear power plant, tliere is always a 
finite probability that a system may fail and lead to an 
undesirable situation. An equipment may fail due to some 
disturbed condition. The transient analysis helps us in 
studying the degree of danger that may occur during such 
failure, and hence gives guidelines for proper instrumenta- 
tion to prevent, such accidents. In case of steam generator, 
large number of possibilities which can lead to an accident 
can be enumerated. There is always a possibility for a 
pipe to be ruptured in the primary or secondary circuit. 

Such a rupture in primary circuit csm reduce the flow rate 
and in some cases lead to a very severe accident called 
loss/ of coolant accident (LOCA). In this situation coolant 

i 

at very high pressure will be exposed to low pressure system. 
This coolant which is radioactive can be hazardous to human 
life. Partial or complete failure of one of the valves 
or pump can also create lot of problems, as it can lead to 
an accident called loss of flow accident (LQPA). All such’*" 
disturbances will effect the pressure and temperature of 
primary fluid. Similarly a reduced flow in the feed water 
pipe will effect the steam output, and hence the power 
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output of whole plant, as it will upset the operation of 
the turhine. In order to ensure that the heat exchanger 
Gan accommodate certain deviation from steady state opera- 
tion and to ensure safe operation if there is an accident, 
the transient analysis of steam generator is a must 

1.2 REVI3W OP PREVIOUS MOUEIS 

In the literature there are no simplified models 
available for simulating steam generator excepting the 
complicated models of the Reactor thermal hydraulic computer 
codes such as R.ELiP-4 [2] ajid PLASH [5] etc. Recently, a 
mathematical model was developed by Arwood and Kerlin [4] 
to simulate- the IT- tube steam generator. This paper used 
linearized equations to predict the behaviour around a 
steady state operating point. The detailed heat transfer 
regimes were also not considered explicitely. Hoeld [5] 
developed a nonlinear transient model for the calculation 
of the dynamic behaviour of a vertical natural circulation 
U-tube Steam Generator (UTSG) . But modern power plants use 
lEUTSG which is an improved version of UTSG and has higher 
heat transfer coefficient. I’into [7] has developed a model 
for steady state simulation of the IBU^SG and Silva [6] 
developed a model for simulating tr ansi aits of an IBUTSG. 

The above two models take various flow regimes into account. 
The second one by Silva [6] follows an imxolicit numerical 
technique for finding the dynamic behaviour of the steam 
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generator. Similar work with some modifications was 
simultaneously done hy two people recently. One was done 
by Silva [9] using implicit method and the other one was 
done by Kar [l] using explicit method. In Kar's work 
there are some numerical instabilities, around 0.4 seconds. 

1.3 PRESENT WORK 

The present work is an extension of the work done 
by Kar [l] for his M.Tech. thesis. 

It consists of developing a model to simulate the 
transient behaviour of an Integral Economizer U-tube Steam 
G-enerator (IBUTSG) utilised in 600 MV/ e Pressurized Water 
Reactor. To give an accurate representation of steam 
generator due considerations have been given to two phase 
flow and heat tra)isfer. 

As fundafiiental equations, the well known conser- 
vation equations for mass, energy and momentum and Fourier's 
heat conduction equation, all of them given in the form of 
partial differential equations have been chosen. Additional 
approximate equations for the thermodynamic properties of 
water and steam, for a two phase correlation, for one and 
two phase friction coefficients and heat conduction equations 
have been established, 

E'or steady state values, the model and computer 
code developed by Pinto [7] have been used. 
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Some important modifications have been made in the 
program developed by Kar [l] , for the study of transients 
of a heat exchanger { IBUTSG-) . The instabilities in the 
scheme have been eliminated. Program was run for various 
possible accidents and some important conclusions have been 
drawn. Varying the amount of perturbation, feasible ranges 
of perturbations, below which the program can be used, have 
been established. By running the program at various time 
steps a suitable (economical) time step has been found. The 
final results were verified with those of [1] and [9]. 

In the second chapter a brief discription of the 
model lEUTSG design and it's functioning have been given. 
Simulation model of the steam genei’ator as well as schematic 
representation of steam generator are also given. 

In the third chapter the mathematical model, 
different functional relationships that were assumed, formu- 
lation of system of equations and it's solution and the method 
of calculation are given. 

In the fourth chapter the important modifications 
that were done for the present work have been presented in 
details. Later, the importance and effect of all the 
modifications are discussed. 

In the fifth chapter the final results have been 
discussed and some important conclusions have been drawn. 
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In Aprpendix I the steady state simiAi^tion done by 
pinto [6] hae been given. 

In Appendix II have been given, brief description 
of various subroutines that were used for transient 


simulation. 
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CHAPTER II 

2,1 DESCRIPTION OP THE SI STM 

A schematic of the Integral Economizer U~Tnbe 
Steam Generator (IBUTSG) is shown in Pig, (2.1). 

This is a natural circulation steam generator 
with integral feed vxater pre-heating which produces dry 
and saturated steam. The steam generator, consistant with 
previous designs, is a vertical shell and.U-tuhe heat 
exchanger. The 48500 sq.foot surface area tube bundle 
provides the boundary between the primary and secondary 
flow paths. Incorporated in the tube bundle are multiple 
tube supports and anti-vibration bars in the cold leg side. 
The entire tube bundle .is enclosed by a wrapper to provide 
controlled distribution of recirculated water. The feed 
water inlet nozzle is installed in the lower shell where 
the steam outlet nozzle is located in the elliptical head. 
The upper internals of the steam generator consists of 
moisture seperator and associated support structure. 

The model lEUTSG steam generator produces dry 
saturated steam by transfer of heat from the reactor 
coolant water to a steam v/ater mixture on the secondary side. 
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2.2 DIFPERBNT GOIPONMTS Ca? THE MODEL lEUTSG 

(A) Primary Side ; 

Tile hot primary fluid coming from the reactor 
enters the lEUTSG at the inlet nozzle and flows into the 
primary inlet plenum. The fluid then enters the upward 
portion (hot leg) of the U~tube, where it begins to lose 
heat throu^ the tube walls of the secondary fluid. The 
primary flow continues into the downward side (cold leg) 
and into the economizer section, 

(B) Tube Metal V/alLls % 

The tube metal walls separate the primary and 
secondary fluids. Heat is transferred from the primary 
to the secondary fluid through the tube vralls at a rate 
determined by the overall heat transfer coefficient and 
temperature difference between primary and secondary cir- 
cuits. The overall heat transfer coefficient is dependent 
upon the tube metal conductance and the film heat transfer 
coefficients at the inner and outer surfaces of the tube. 

(C) Secondary Side ; 

Just above the hot t cm plenum there is an economizer. 
This is nothing but a closed container covering portion of 
the cold legs and seperate them from the rest of the second- 
ary side, Peed water enters the integral economizer section, 
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where baffles direct the flow upward across the tubes 
(counter flow relative to the primary liquid in the cold, 
leg). This brings the cold feed water into a region where 
heat is transferred with the cooler primary liquid. The 
baffles in this section promote turbulance and therefore 
increase the film heat transfer coefficient of the secondary 
liquid. 

The feed water mixes with recirculating water to 
form a slightly subcooled mixture that is forced upward 
around the tubes of the hot leg section. The temperature 
of this mixture increases as it takes heat from the pri- 
mary liquid until the boiling point is reached and steam 
is formed. The steam water mixture continues flowing upward 
part of the economizer section. The steam water mixture 
enters the seperators where the saturated liquid is seperated 
from the steam and returns to the dovmcomer while the steam 
is collected in the upper part of the steam generator shell. 

The IBUTSG design is a modification of the U~tvibe 
recirculation type steam generator (UTSG) used in many 
operating PWR. In order to increase the heat transfer 
effectiveness of the system a change in the design was made. 

The IBUTSG is same as standard UTSG excepting -tiiat in the 
UTSG the whole feed water is introduced into the downcomer 
and there is no preheating (economizer) section. As the 
recirculation ratio is very high, the circulated water that 
enters the heated section from downcomer is near the saturation 
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temperature. Hence there is very small temperature 
difference between the primary liquid in cold legs of the 
tubes and the secondary water, in a UTSG. This will lead' 
to a rery poor heat transfer rate at the cold leg section. 

By introducing colder feed water directly in this region 
the lEUTSG design reduces the temperature of secondary 
fluid and hence increases the temperature difference. 

This in turn, will increase the heat transfer from primary 
to secondary liquid. 

2,3 MODEL OF THE STEM GENERATOR [6] 

The steam generator model used in the present 
study is shown in Fig. (2,1). The steam generator is 
represented in two sections, preheater (economizer) and an 
evaporator. In reality, the steam generators of the PWR 
reactors have U-tube geometry. The present model simplifies 
it to parallel tube geometry but retains all the significant 
factors such as the total heat transfer area, total amount 
of heat transferred, preheater etc. Lack of detailed data 
about the internals of the steam generator such as the 
pressure drop due to baffles, quality of steam at the 
outlet of the generator before the steam seperates etc. is 
circumvented by the pressure bad. an ce and heat balance 
equations as described in the later sections. 
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As indicated in fig. (2,2) the coolant from the 
reactor outlet enters the primary side of the generator as 
a hot stream in the evaporator section and leaves the 
generator from the preheater section as a cold stream to 
re-enter the reactor. The heatflow from the primary side 
is conducted through the tube walls of the steam generator 
and enters the secondary stream. The figure also indicates 
that the secondary stream leaves the evaporator section of 
generator as vapor and part of the steam is recirculated. 
Some make up feed water enters the preheater and mixes with 
the Secondary stream after the preheater. 
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CHAPTER 3 

3.1 mathematical MODEL; 

1 : simplified mathematical model and computer code were 
developed hy A.O, PINTO [8] for the steady state analysis of a 
steam generator utilized in a 600 Mwe pressurized water reactor 
nuclear power pl/ints. As a reference to present work the reprint 
of the above paper is included as Appendix 1. 

3.1.1 Eundamental Conservation equations ; 

Some assumptions are made to derive the equations that 
are given in this section, 

1. The flow is vertical 

2. There is no transfer of work across each control volume 
indicated in fig. (2,2) 

3. The area of flow cross section is constant 

4. The effect of gravity is negligible. 

Thus the three fundamental equations of mass, energy and 
momentum conservation become 

Mass : ^ ^ ” If (3.l) 

Where W = pAV 

8 (pe) • 3Aj 3 (We) 3 ( PAY) ^ " . 
Energy ; A = q (3.2) 

Where e = u + V /2 + gZ is the energy per 
unit mass of the fluid. 
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Moment'um i 


at 




- A pg 


a 

a? 


(v/v) 


( 3 . 3 ) 


Because of the 4^^ assumption the terra gZ in the 
expression of 'e' and A g in equation (5*3) can be 
neglected. 


3.1.2 Pinal systaii of equation for tranaient simulation 


at az “ a z “ 9Z 


A 


lizzl 

a t 



a 


(Wu) - PA 


+ V 



(3.4) 


(3.5) 


A 


at 


ay/ • 

a z 


(3.6) 


3*1.3 Integration of conservation equations : 

Integrating the above equations over a control volume i 
hounded by junctions i -and 1+1 with respect to Z lead to the 
following equations s 


(a) Ihe equation of morn en turn cons er- vati on gives: 

, i+1 i+1 i+1 

ap, 


i+1 


i_ / , a. = - / / ^a.-A f ^ 

X X X X 


Let us define 


/ PA V da 


i+1 


dZ 


1 


dZ 
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With this definition we get 
dW. 

4 Z = -A(P - P^) - [W - W . ] 


P 

Ki 


But W = PAV, i.e. , W = V/ /PA, and 

- 2 


W, 


^Ei 


f . 

3 . 


AZ 


2 AD 


e 1 


So the momentum equation takes the forms 

1 * _i_ (!i _ 

/ ' l'i A j-ac ' 


«i = h 


A. A Z P 


) “ f. 


2 

^i 


'i+1 


2ADe Pi 


(h) The equation of energy conservation will give; 


(3.7) 


i+1 


i+1 


i+1 


i+1 


9 

TT / Apu dZ = 

t i i “ i - i i 


i+1 


I 5 ^ dz - J PA^az + J v'j^z 


We now define 

i+1 

= / A Pu dZ 


and 


Qi 


= ,i q az 

X 


i+1 rsTT 

! PA ^ dZ = P. A(V. = P. 

i 


and 


i+1 


9P 


K 


I V -y| az = v^Pjj^. 


'1 


3 


^ 2a2 


With these simplifying assumj)tions ve get 


- ■«H-l 




w. 


0+f. 


A Z 


«1 


i+1 


^2A^ B^P. ^ 

t/ X 


(c) The equation of m(a.ss consei*vation gives; 
i+1 i+1 

it / oA dz = - / H a^ 

X J- 

Let us define 


A Z 

(3.8) 
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i+1 

I A P dz 

P “* ^ ^ 

i i+1 

/ A dz 
i 

With this definition we get 


* 



•’1 

A. A Z 





(2.9) 


Now we finally rev/riting the above equations 


•i. 


rkz P'l - 


(3.10) 


A 


Wl 


w 


Q, + - Pi C22( 


i+1 


i+1 


W. . 

-i) (5.11) 

Pi 


+ 0 


A 2 

23 2A^ Dp ^ 

K* -L 


^91 ^ 

-‘^+— fp. — P. ) ^ 

A 2 '‘•^i ^i+1^ 


r, 2 


h ®i 


1 ^ JL+J- \ 

2AD.P. ‘’1+1 


^'1 ‘^i+l'l+l 


e 1 


(3.12) 

3,2 Relationship between Junction Variables and Control 

Volume Variables and Corresponding Derivatives Among ©lem; 




1 

2 


(V/i+Wi+i), i,e,/ >^i+l = 2¥i “ ^i 


p; 

1 


1 

= 

( + A+i)> 

i. e# , 

A+i = 

2®! - '‘i 

11 

'( ”1 ‘'i+i^ ’ 

X « G 9 9 

*^i+l ” 

2 - "i 


D ef ini tl on of sp eci fi c eint ernal en er gy 

Apu dZ 


i+1 
/ 




U. 


i+1 

/ Ad’dZ 


A. A Z p. 


(3.15) 


i 
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Some of the derivatives obtained from the above relations 


9 u. 

2 

M“ji 


A. A Z P 


u 


1 


: 


a p . 


A. A Z "p 




p ’ 


?P 


U’l ^ 


a p 


a p 

i+l _ 
“a P. 


- 1 


5 u 


1+1 


= 2 


9u, 


9 U. - 
1+1 _ 

9u. 


3 ¥. 


i+1 _ 


9 \h 


'■iTq~ 


5,3 Functional Relations that are assumed 


fi = 


Pi = ^(^it ^i^ 


^i+1 ^^^S'.+l' ^i+1^ 
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^i+1 

= f(%, 

\) 

^i 

= 

\) 

*^i+l 

= f(Oi. 



= 

Pf) 

•H 

= f(Uj,, 

Pi) 


3.4 SOLUTION OP STSTEW OP EQUATIONS; 

The v;hole length of the channel is devided into 
different control volumes. Tlie three conservation equations 
consists of three unknowns, porching [l2] has given a scheme 
which can be utilised for solving hydraulic networks. Similar 
scheme was derived in ref. [9]. Porscliing's scheme is veiy 
powerful for nonlinear problems. This can be used explicitely 
and implicitely as v/ell. The explicit scheme is going to be 
used in this work. 

3o4-,l Scheme suggested by porsching; 

Por a given initial value problem of the type 

f = P(t, y), y(o) = y° (3.14) 

One can use an explicit method given below 

[l - h dP (t^, y^)] A y^"^^ = hP (t^, y^) (3.15) 

^ yn4*l ^ y 114*1 ^ 

3.4.2 Scheme derived in ref. [9] 
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In the matrix form it will he as follov/s 




a 


5.4.3 Derivation of present sclierae used; 

[l S At] At 

^ n ■“ 

The given differential equation can be written as 


li = £ (*■ y) 

.. . . n +1 _ -^n+l I 

i.e, , Ay = P At 

= . P^ + - P^j A t 

= , ( + A At 

If we approximate 

dP _ A jjJ^+l 

^ n 'Ay’^+l 

then A p ^+1 ^ n ^ ^ 

So we now have 

Ay^**'^ = [ P^ + -S ^ A t 

- ay - 

i,e., [I ~ -p A t] A = P^ At 

UJ — 

The matrix representation of this equation will be 


( 3 . 19 ) 
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• 

a 

1- At V 
^ Pi 

# 

8 p • 

_ At -— 

^pti. 

1 

-At 

• 

tfl 

a w^ 


AP. 


* 

^i 

^•&i 

- At 

a\ 

30. 

1- St -i 

aS, 

"• At 

aWi 


Atii . 

; = 

. 

4 

aw. 

- At ~. 

4 

a ¥. 

„ At ^ 

■ au. 

1~ At 

A 

aw^ 


AWi 


Wi 





1 

— 


— — 


A t 


A X = B 
B = P . At 


( 3 . 20 ) 


5»3 Description of various matrices of equation ( 3 * 20 )' 


3 . 5.1 Matrix ' P ' : 


= aTTz - Ah-P 

- '''l ..1 Wh 

^2 = ( Ha ,! - “22 


f - 


L W? A Z 

- I ' C 2 'z : r ^ ) 

A "^ D ^ p .^ 

Q J - 

( P .| - P ..-,) ^ 


^^ i + l ^ i +1 


• i + 1 - 


2 -^ Vi ^ 


3 . 5.2 Derivatives of ? ¥ j _ 


a p . 


a p . 
1 


= 0.0 



23 


Hi 

9 U . 


= 0.0 


9 P , 




= -2/A, A 


Hi 

a Pi 


H 


A,i 


Hi. 

9 P 4 


2 ^^-i ,1 

+ Cgo Pi ( 7 ^^ 

22 1 Pi^.i 


au 

3 U . 


i 




Ha.-i (--r 


V/. ,T 
^ 1+1 

^ A 

w. 

_1.) + 2 

0 • 

"i 

'■'i+l T 

^i+1 

* 1 

'’i 

0 

ro 

VJ^ 

H-bl 

iH‘ 

Wi^ A z 





• au, 

. ”) 
9Ui 

- %+i 

■,0, • 


au . 
3 ij . 


’*‘i Ali f ! i .+ X . - !ii 
°22 au^ ' atJ^ '’ 1+1 “l 


30 . 


sti 


9 VL 


9 p . 
1 




Hi±i 

9 ?L 


G 




3 VJ . 


i+1 


^ i+l 


9 VL 


o „ I , 3 w 2 a z 

+ zr^ 

2 a 2 Pi 
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A ap ,- 4.-, , 2 ^'^ i+lHiJi 

/■ „ ,.; k ± P ) + -i — 1 — « + 2 

^ “ 2 A D„ -A. P 4 .1 


a ?. 

1 


e 1 


i-M 


l!i 

3 Wl 


'* '’^ i+l ’■^ 1+1 


:!iA 

1 P. A /s Z ^i^.1 
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3.5.3 


3W. 

1 

StJ. 

1 


"21 


A2 


S U. 
1 


Derivatives of pressure And Temperature s With, respect to 
Independent Variables Density ( f ) and Internal Energy (u): 

(i) Derivatives of pressure; 


3 p . 


1 r 14-1 

2 *• . - 

9 p . 


O 

+ ] 


- ap. 

1 [C-^) 


2 ■- ^ 9U 


3u. ^ au. 

(-^)„ C-“)i 


, P \ ~ .'ll, ''u. 

1+1 i+1 a U . ""i 9 ""i 


9P. _ 

+ Ct+^) 


3p. . 


^ \+i ^i+1 a i: 


- 9p. 

^ [ 2 ( 




1+1 


- 2 (^Op 

^\+l 1+1 i 


3?. 

(“) 
a Pi 


8P. 

1 


U. 9 P. 

1 ( 1+1 \ . 

^^i+l'^^i+1 Pj^ \+l i+1 


9P. ,, 


./ 


/O 


, 9P. , 

1 [(_^) 

^ 9U> /a^. 


jO 


ap. ^ 
1+1' 


(3.21) 


2 Vl ^1+1 Av, "^i+l i+1 ^xs. i 




T .1 

^^i+l i+1 9Uj_ 


au. 3u. 

(_yi) (^)] 


su 


3P.. 


= %[(^) 


";p 

*^^+1 i+1 


A A Z Pj_ 


2 
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9P, 


au, 


A A 


p. 

1 


aP. - 
(-i+i) 

*"1+1 


(ii) Derivatives of Temperature: 


( 3 . 22 ) 


hSi 

II 

f(Ui, \) 



a!'. 

(~tr) 

^pi 

3T. 

= (“)5. + 

D i 

32. 

(*^)p. 

3Uj^ 

3U. 

(--) 

aPi 


aPi 1 

3 2. 

(“r^)T5, 

3U. 1 

U'. 

(- zr) 
Pi 


32. 

= .. 

3 2. 

=■"1 ^ 

u. 

(z^) 

Pi 




3T. 
!• 




(^) = (-^-0 
au. 9u. ^i 3U. 

1 11 


aT. 

= (-^h. 


1 


9u^ 'i A AZ.P^ 


(3.23) 


(5.24) 


3.54 Elements of Matrix 'A' 

a._ = 0.0 


12 


^13 


^ ^A. A 


a 
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80i 1 Pi SU, 'iJ Pi 


\ 


- 2 Coo '~r^s + 0, 


'22 
W-- 




7T 


+ (|^>u - J ] 

22 oji+l %+l ii "^^i+l ^i+1- 
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= 1 


sT. 2V/. . 

^ 1+1 


k & Z 7 i AZ S . 


. 2 Opp , Cp, f. W.® AZ 

A m r ^ n t <.<-1 ^ 25 1 1 H 

^ 1+1 2 “2 ^ ^ 

* i+1 P.. 

C ? -L 


/t / _ 1+1 1 \ 1 

0 ^ 

i+1 Ji AAZS, 
„ i 

^i+1 ’Sl+lt 


^ 2 Cpn A • 3p. u. 9 P,^t 

A T [ { (r~)^ + (ni^')p > 


w 2 

2 . A . D ^ ' P .‘ 

v5 1.1^ 


^ ^^ 1+1 ^ i +1 

A - <’!+!• «" 


r 2 Gp-, 3P. T 

^ T [ (f5^) 3 

( aZ )". ^^ i +1 


1 . at r ^ i-A ^ i+ll 

A^,,AZ 


VBOffOR 'X': 


x(l) 

- aV 

= 

^ T 

•«» P* 

X 

X(2) 

= M-j/ 



X(3) 

- 'f 

= A 


- T 
- W- 


3.6 Method of Oalculationt 


Steady state values have been taken as initial values 


for the given problem of type 

y = P(t, y), y(o) = y^ 

Introduoe a step change in one of the given 
variables : u, p, Vf, P at time t = 0 and at 
the inlet. 
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Solve for the increments at any 


time t equation (3* 20 ) at all control volumes explicitly. 

Then add these increments to values at previous time 


step to get the 

values 

at new time step 

7 T+1 

- T 

, A - T 

Pi 

= Pi 

+ A Pi 

fj. 

= U. 

+ Au.^ 


1 

1 


- T 

= 

+ A-Wi^ 


Knowing the density find cut the specific volume at 
the new time step. 

Finding T, H, X ; 

Oorapressed water i Here isothermal flow is assumed as the 
temperature change from one time step to another is negligibly 


AnT-hI 


small 


i * e . ; 

(4-)" 


i>i 

So 

5. 1+1 ^ 
1 

Then 

T+1 

H. = 

1 


^InT' 


U,. 


^i 


T+1 

i 

T+1 


p. 

1 


. C 


22 


% 


T+1 


Using the table of temperature as a function of pressure 

. rn ijn 

for compressed Water temperature Tj^ can be calculated. 
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Saturated water ; Pressure, quality can be calculated as a 
function of internal energy and density 

X = X(P) 

P = pCS^, u.) 

Knowing pressure corresponding saturation temperature can be 
calibrated from the table, 

Then enthalpy is given as 

H = % + X(H^ - ) 

saturation enthalpy of liquid 
and vapor corresponding to the above pressure P, 

Definition of some constants used in this chapter, 

2 • 2 

Got = X 5600^ X 144 = 0.600444 x 10^ 

° Hr^ ft^ 

O22 =144 ^ X 0.001285 = 0.18504 

2 

C,, = — X 5 x 0.001285 MH - y .. - 0.308172x10"“ 

25 So 360 o 2 Sec^ 

5,7 Flow Chart i 

The flow charts developed for steady state simulation 
and transient simulation are given below. 




ARRANGE VALUES IN TABULAR FORM 


DETERMINE THE AXIAL DISTRIBUTION 
OF TEMPERATURE AND QUALITY OF 
VAPOUR IN THE EVAPOITATOR 


GEOMETRIC A! 


CALCULATE THE POWER GENERATED 
m THE STEAM GENERATOR AND 
COMPARE WITH INPUT POWER 
RECEIVED FROM REACTOR 


OF PRESSURE % 


EXCEEDED MAX. NUMBER 0F;::>- 
INTERATIONL--^ 


CACULATE THE PARAMETERS OF THE 
PREHEATERS HEIGHT^ HEAT FtUX i 

pressure^temperature etc.. 
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CALCULATE P-Ur/w AND U Qy IN Al l, 
COMTF?Ol VOLUMri:S AMD AT ALt 
JUNCTIONS AT t 0 







INITIAir^E THE CONTROL VOLOM'c IN 
WHICH ONE-OR TWO PHASE FLOW 
STARTS AND END 


v? 

r^S'«.T^i’'-Tr4 .it-n, ^ ■ ! k g. ■ r r « r. v* 7 » « v* .(<•*.- >. |i * . , i ^ • r r-- >'-•>’ '-tb • 1 7 - /C' ?V' » •* 

INITIALIZE THE T!?v)E AMO CAClJI.ATf:' 
THE TIME STEP 


DO 10 I .. 1,MIT I 

Urr^TOTAL NO OP STEPS | 




'V,. ' 


UETERMiNE WHICH CONTROL VOLUME HAS 
ONE- OR TWO-PHASE FLOW AT 11MH >!>’ 


Tall' a m6¥o ' aTFd ” 
^HIAr,.F„LUX 


DETERMINE THE AMOUNT OF HEAT TRASPERR* 
IMG IN EACH CONTROL VOLUME 


DETERMINE THE AMOUHT OF HEAT 
TRANSFERRING IN THE PIIEHHATER 

7 In I fc p > »pipi r i T B t i^ t | M- ii rT-T-rtri^ nrT-r It • r --— ^ - -r- -^ - — ■>— 

i 

I 

. I II I ^ I III I ■' - - - m.H. '"'**3 

DETERMINE THE MODE OF FLOW IN THE SECONDARY CIRCUIT J 



DERIV 

CAlCULATE(0T/3f \j OT/3U)^ IN EACH OF THE CONTROL VOLUMES 
OF THE PRIMARY OF EVAPORATOR -' ' ' ^ '■ ■' ■ 








J JCrK. ^ 


WjlTlilV® 

9AMraiA'j:K{;i'i'/ejS)vjc (QfAmh xw m\on m im GO]f?K,ioL v()X.MMm3 
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GH^TER W 

4.1 M OD IP IC ANIONS PIADB IN THE STEADY STATE SIMlHiATICN 

The steady state simulation developed by pinto [6] 
has some limitations. There the only thing that was fixed 
in advance was total heat flux. The way in which simulation 
has been done is given in Appendix I. An assumption which 
was made in it is that the steam flow rate (Wg) Is equal to 
quality (X) times the total flow rate (W^) at the exit and 
is always fixed. 

'^S = *exit • ''Total = constant 

So in order to match the total heat flux one can 
have many alternatit'es of total flow rate and exit quality. 
There can be a higher quality mixture vrLth less flow rate 
and vice-versa. But this particular program is converging 
for only initial approximation for quality which is around 
0.4. It also always converges towards the same quality which 
is 0.39. 

As mentioned in ref. [ll] a crude measure of adequacy 
of secondary coolant flow through the tube bundle is the 
"circulation ratio" which is defined as weight of viater 
circulating in the exchanger per unit-v/eight of steam 
geherated, A small circulation ratio indicates a large 
fracti on of steam in the tube bundle and the possibility of 
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poor heat transfer on the shell side. In addition to this 
a low circvtlation ratio leads to flow instaMlity. In 

the past, anits had been designed with circulation ratios 
such tljat the steara quality at the exit of the bundles did 
not exceed 20 to 25 '/» . 

A low circulation ratio can lead to flow oscillations 
that cause periodic uncovering of a portion of the heat 
exchanger surface. This ‘chugging* phenomenon reduces the 
heat transfer efficiency, and flow oscillations of sufficient 
amplitude can lead to large fluctuations in the water level 
or steam flow rate. 

Experience shows that flow stability can be achieved 
by keeping the steixfri quality in the bundle low, or equivalently 
by keeping the recirculation ratio high. 

So it was tried to reduce the quality at the exit 
by assuming low qualities around 0.25 to 0 , 3 . But the scheme 
did not converge. Many modifications have been tried but 
they failed. At last the total flow rate \ma made constant 
and the steam flow rate and quality left open to converge to 
their own values. By this a control over total flow rate was 
achieved. As expected this has worked well and quality 
started decreasing with the increase in flow rate and vice- 
versa. Einally a quality of around 0,3 has been chosen for 
transient simulation, ; 
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4.2 MODIFICATIONS MADE IN THE TRANSIENT SIMULATION 

In earlier formilLation the results exhibit some 
instability which produces large variation in the parameters 
and almost all the results are converging to abnormal limits. 

As a first step the whole problan was reformulated 
and compared with previous formulation [l]. It was found, 
there is no error as far as the mathematical formulation 
is concerned. Though it was mentioned that the instability 
is purely a numerical instability, no attempts were made to 
pin point the problem. Some of the possibilities are ; 

1. There may be a big calculation mistake in the 
program. 

2, The scheme used may not be logically correct. 

5. There may be typographic errors in the program. 

If everything is alright as mentioned in the previous 
work, it is purely a numerical instability, and a better 
time step will automatically eliminate the instability. 

As it was mentioned that the instability could be 
due to large time step, the first trial that was made with 
the program, was to study the program vrith reduced size of 
time steps. When it was Ijried with time steps less than 
previous time step (0. 001} , instead of converging to bett er 
solution it s tar t ed diverging. At this stage if one is very: 
0 onf ident about the oorr ectness Of other things - it can b e 
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concluded that this method will not converge for the 
present problem. 

In order to see the effect of the scheme tvfo simple 
schemes were tried. One is Euler's method (most simple of 
all) and the second one is Fourth Order Runga Kutta Method 
(very accurate of one step methods). 

The Euler's method started blowing up immediately 
after four to five iterations. It has given an idea that 
the accuracy of the method is also very important for the 
present problem. It was confirmed when fourth order Runga 
Kutta method run smoothly. But the instability persisted 
there. This has revealed the truth that there must be some 
calculation or logical mistake. 

Proof reading of the program was done several times 
as a result of which some typographic errors were found. 

At two-three places the brackets were placed at wrong 
location. But the mistake is such that the computer can 
not find out the error during compilation. Ofcourse this 
correction did not change the solution very much. 

The most crucial calculation of the program is the 
calculation of increment function from equation (5,20):* 

This involves the inversion of a 3X5 matrix at all the 
control volumes ahd at all the time steps. It is very well 
known that the matrix inversion is so sensitive that it 
reacts greatly even for smaH fluctuation in the yalues of 
it 's elements. 
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So the accuracy of the elements of matrix A (3*20) 
must be ascertained* As it involves several derivatives of 
pressure and temperature, accurate evaluation of derivatives 
is needed because of the fact that the calculation of the 
derivatives is only possible through the interpolation of 
curves given in steam and v/ater tables [ 13 ]* Ojn verification 
it was found that the interpolated values of the derivatives 
were very erroneous. The new values of the derivatives have 
been interpolated and introduced in a new subroutine (DERIV) 
for finding derivatives, in the place of the subroutines 
DPRBU and DPRER. 

Not only the value of derivatives but also the method 
of expressing (finding) partial derivatives from functional 
relations is important. So the partial derivatives have been 
expressed as sho-sm in section 3*4. 

Here is an important modification which has. almost 
nullified the instability. In the previous work there is an 

JU t 

assumption made while finding the derivatives of U with respect 
to p^, u^, The term Q in the energy equation (3«I1) is 

assumed constant with respect to average density ^ and 
average internal energy u^. The simple expansion of Q explains 
why this assumption is wrong, 

^ “ ^A,i ^^0,1 - 
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^A,i all heat transfer coefficient between primary 

and secondary at control volume i. 

Average temperature of the fluid in control volume i. 

Here the average temperature of secondary fluid is 
a strong function of internal energy and to some extent 
density, 

-3U. tu. 

The elements — of matrix A ( eqn. 3*20) 

3 T. 3T. 

involve the terms(H. . • ■— — and ( H. . • — which are 

8(5. 3u. 

1 1 

not negligible. Hence it can be concluded that ■ ' ■ , are 

3 P . 8 u. 

p1> 

not only non-zero but also significant. 

With the above modifications the instability has been 
reduced very much. Yet there is no convergence at all as far 
as mass flow rate is concerned. Except the mass flow rate 
all the other variables are converging. But the quality 
rises to a very high value around 0.7 to 0.8 and converges 
at that point. Whereas the steady state quality is around 0.4 
Just by inoreasing the energy content of primary liquid by 
5 /. it is highly impossible physically to get a quality jump 
frcmi 0,39 to 0.6^ increment of about lOG/, ' for a small 
varihtibh in the primary side. At the same time the flow rat e 
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rate is hi^ly fluctuating, This has given a clear indica- 
tion that neither the mass nor the energy is conserving 
properly. The previous work did not "bother to verify the 
total heat transfer rate with the time, "When it was tried 
to get the printout of heat transfer rate, surprisingly it 
was reducing to very low value with the time, which was as 
low as 50 % less than the steady state value. 

At this stage it was tried to verify the logic very 
carefully. During this verification some logical errors 
were found, 

THE LOGIC USED IN THE PREVIOUS WORK ISs 

After deviding the channel into N number 
of equal parts, the conservation equations have 
been solved f or the increments in the first and 
Nth control volumes only and it was assumed 
that the average of these two increments can be 
taken as total increment from bottom to top. 

Also this total increment can be distributed equally in all 
the control volumes. So in the previous work they did not 
tried' to include the effect of the values of the intermediate 
control volumes. Also it was assumed that the distribution 
of various parameters along the length of the channel with 
time is parallel to the initial hi st:^ibutiGh (steady state 
di S tr i button) ■ i 
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It seems that these assumptions are wrong. The first 

assumption is inGOrr<^ct, because it involves gross error, and 
it has to he corrected. Though it is not the correct approach? 

there is no harm in considering the second assumption. The 

increments here are of the order of less than 0.1 % of the 

actual value. So it is reasonable if the incraaents, which 

were found out at different control volumes explicitly, were 

averaged over the length of the channel. In other words the 

increment at any time step is constant and is equal to average 

of all the individual increments. 

Mathematically if is a small increment of a 

function at any control volume i, then the logic of previous 
work is as follows ; 

-T 

Pind A \ ' and A 


-T 


Inoremont in the first control volume at time T, 


-T 


- Increment in the Nth control volume at time T. 

-T ~T 
-T 

■Qj. _ 


2N 


At any control volume i 


^T+1 


-T 

Y: + DY 
1 
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The correct logic that should be usedt 


Find Ay^ 


i = 1,2, (explicitly) 


-T -T 
Y -|- A Y 


Another logic that can be used reasoneibly is 


Find AY^ 


i “•'1,2,*4 *«,jl 1 ( exp Ic. c x tiy ) 


I a: 


-T+l 

Yi 


-T 

Yj_ + FY 


There is a semi-eXplicit scheme v^hich can be derived and 
is very powerful for some problems. 


Find A 


Using the inlet vlluos of control volume 
i at time T+1 and exit values of control 


volume i at time T. 


Find Y 


-T “T 

y 4 - A Y 

^i i 


Using the expression 


„T+1 , . yT+l 
* !i±i 
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yT+l _ “T+1 ..T+1 

M+l “ 

Knowing the values of junction 1 at time T+1 one can 

go forvrard to get tho exit values. 

After this modification, another logic step ims 
found which was uni’ eas enable. There are some statements 
which state that whenever one of the variable become negative, 
after adding the increment, subtract back that increment and 
proceed with the previous value. This is definitely a suppre- 
ssion of error vdiich is existing at that point. In the basic 
progrt®, for which some results have been given, the removal 
of these statements produce results which are not correct. 

In other words, tho program is not at all running properly. 

While searching for the reason for it, it was found that flovr 
rate becomes negative after some time. So the program is 
becoming unstable. This is definitely an indication of error 
and one has to stop at this stage to search for the error, 
because the increment should always be within the limits. 

After the present modification this type of error is not 
reoccuring. 

Coming back to the heat balance, simply speaking, 
the heat transfer rate is approximately equal to = V/ ill. 

Where W is the average flow rate in second cry of evaporator 
A H is Enthalpy change in secondary of ovap orator , 



On the other hand flovr rate can be written as 


W = Q/ 6 H 

where, Q is the heat transfer rate calculated from the 
correlation of various heat transfer regimes along the length 
of the channel. As the flow rate is highly fluctuating it 
was expected it is mainly due to erronious calculation of AH, 

As shown in fig. (2.2) at the inlet of the secondary 
of evaporator a jet of recirculated water at certain conditions 
will mix with another jet of feed water from the preheater of 
the secondary side. To start with the calculation of var'ious 
parameters along the length of the channel one has to know 
the inlet conditions. For the secondary of evaporator inlet 
conditions are nothing but the conditions of the mixture. 

Because of the lack of details about the interAe3t.ls of the heat 
exchanger, it is difficult to make a momentum balance or heat 
balance for finding the conditions of mixture at the inlet so it 
is necessary to assume appic'oximate values for the parameters. 

It was assumed that the total enthalpy and specific volume of 
the mixture consists of their individual values in proportion 
to their mass flow rate, i« o» 
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^2 ~ flow rates of two jets 

Hg ~ Enthalpies of two jets respectively. 


This assumption is reasonable. Before coming to this 
point we need to state that there is another mistake committed 
in the previous work while calculating liquid, enthalpy of 
recirculated water. Here is the formula that was used 




^ exit ^ 


H 


exit 


If one thorou^ly goes through this expression he 
can be convinced it is not correct. 

Expanding the above expression 

(1 - ^exit^ ^^llq “ “ ^exit^ ^exit 


Vj^hen quality at exit i. e. ^gxit ^ then ^ which 

has not got any meaning at all, because when quality reaches 
one, recirculation flow reaches zero. So, will become 

zero but not the enthalpy content of liquid at that point. 

Also if = 0.5 /than = « which is / 

improper so if the qualiiy at the exit reaches very close 
to 0. 5 the enthalpy content of the redirculated water reaches 
very very high value, and ultimately at 0.5 bsoomes infinite. 
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It can be reasonably assumed that the mixture 
is in thermodynamic equilibrium. As the mixture is at 
saturation state, the saturation liquid enthalpy can be 
taken as enthalpy of recirculated water. V/ith this 
modification, proper heat blance has been achieved.. But 
there is no convergence as yet. Plow rate is still 
fluctuating. 


In the previous work no assumption of perfect 
control of feed water was made, i.e. The feed water rate 
is assumed to be constant throu'^ out the transient. The 
disadvantage of this assumption is as followsj For instance, 
if all of a sudden the energy content of the primary liquid 
is increased and kept constant, then, the heat transfer 
rate from primary to secondary circuit will increase. So 
the quality at the exit of secondary circuit will be increased 
The increase in quality will increase the steam flow rate and 
reduce the recirculation rate. As the total flow rate in the 
evaporatory is the sum of the recirculation and feed water 
flow rates this total flow rate is going to be reduced. 
Because the feed water rate is still same as steady state 
value. AS the less amount of liquid recieves the same higher 
amount of heat flux, quality will still increase, this inturn 
will improve the steam flow rate further and decreases the 
total flow rate in the evaporator. If this process continues 


ther e may 
become 1, 


be a possibility that the quality at the exit may 


and hen ce 


then 0 will be no reciroulation at all. 
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The feed water that is sypplying will be directly converted 
to steam. This may lead to burnout condition also. In order 
to avoid this there should be some control over feed water 
rate. A possible elimination of this sort of trouble is by 
feeding more feed water. That is if it is assumed that 

^feed “■ ^exit ^total ^steam 

This means that the feed water rate should be adjusted such 
that the feed water rate can always match with the steam 
flow rate. 

One more mistake was found at this stage. It was 
another wrong assumption made in the previous work, that the 
pressure at the inlet of the secondary of evaporator is 
equal to pressure at the exit of the evaporator at the 
previous time step, i.e. PSJ(l)^^^ = PSJ(EXIT)^» As there 
will be a pressure change of only about less than 1 PSI 
between any two time steps, it is improper to assume this. 

This means the pressure drop frcto inlet to exit Is negligibly 
small. This will obviously affect the flow rate considerably. 
So it was tried by assuming the pressure at the exit of the 
pp@-.h 0 ater equal to pressure at th© entrance of evaporator. 
This way it is giving better results. But still not very 
oonvinoing. It was assumed that the pressure should be 
inbetweeh these two. By taking the average of th ©se two 
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pressures still better results have been obtained, 

Then the program was tested at varying time steps 
and for varying step changes. The effect of these variations 
have been discussed in the last chapter. 
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GHiPTBR T 

5.1 RESULTS AND DISCUSHTOW 

This chapter consists of the discussion of the results 
obtained after making all the modifications given in chapter IV. 

5.1.1 Steady State Values ? 

The steady state operating values which are used 
for transient simulation are given in. Appendix III. In the 
computer code developed by Rinto [7] a lower total flow rate 
and a higher quality (0.4) have been used in the secondary side. 
But in the present case some what higher flow rate and lower 
quality (0.5) is used. By reducing the exit quality on 
secondary side, the region in which film boiling exists, is 
reduced. This is because of the fact that the film boiling 
is unstable compared to nucleat boiling. Besides this, the 
heat flow correlations are not very accurate in film boiling 
compared to nucleate boiling# Plow is also more stable at 
lesser qualities compared to higher qualities. 

5.1*2 Transient Values. ; • 

The results of the transient analysis are presented 
in graphical f orms from fig. 5* 1(a) to fig. : 5.4(f)* 

(1) 5 *4 inoremeht in internal energy at the i^ 

priDflary circuit: ! The figures 5*l(a) to fig* 5.1(f) 

. ' o bhsi s t s of thi s c as e . ; 
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With, the increase in internal energy on primary side 
there will be more energy available at any time in the primary 
liquid. This will increase the heat transfer rate from 
prjjaary to secondary liq[nid» As it can not transfer the entire 
internal energy to the secondary side the internal energy at 
the exit of the primary side will increase. The temperature 
of the primary lic^uid at the exit will increase, as some of 
the extra energy avsd.lable in the primary liquid can be 
utilized in increasing the temperature of the primary liquid 
itself. Coming to the secondary side the extra heat supplied 
can be utilized in two ways 

1, It can increase vapour content of the secondary fluid 
and thereby increase the quality and flow rate. Here 
again, there oan be two types of explanations for 
pressure variation. One, with the increase in the 
quality, the corresponding saturation pressure may 
reduce and hence the temperature. The second one, with 
the increase in flow rate the pressure drop has to be 
increased and hence the eaclt pressure will reduce. 

As the mixture is at saturation state this will 
reduce the corresponding temperature at the exit. 

2, The extra heat can be utilized ±n heating the mixture 
to a high temperature. With this the o or responding 
saturation pressure will also increase. Here unless 
some extra vapour forms the qua;iity will Ppme 
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(il) 5 / increase in mass flow rate at the inlet of the 
primary circuit (figures 5.2(a) to fig. 5.2(h)) will 
represent this situation. 


As the flow rate was increased, the heat transfer 
coefficient from the primary fluid to wall increased, "because 
Reynold's number increases with the increase in velocity. So 
the heat transfer rate from the primary to secondary liquid 
will increase. 


On the primary side, in order to accommodate extra 
flow rate the pressure at the exit will decrease. As it can 
not transfer all the extra heat due to extra flow the temperature 
at the exit will increase. The flow rate at the exit is sli^itly 
increasing to above 5 % level and finally stabilising at 
5 y. level in a short period. 


On the secondary side, the Increased heat transfer rate 
can produce two effects as discussed in the previous case. 

Here in this case, the internal energy and quality are increasing 
and the pressure and temperature are decreasing. This means 
most of the energy is utilized in improving vapour content 
or quality. 


Mass flow rate W ~ S ^ ^ ^ 


av 


With' the increase in quality, average density will decrease as 
more vanbur forms. At the same time average velocity will ; 
increase as vapour v el ooity is more compai^^ft»1iO, iiaui4'.*%*&Gjifthe 




AllffSi 
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mass flow rate will effected according to the dominance of one 
of these. If the variation in density is less predominant then 
the mass flow rate will increase and if the density variation 
is predominant mass flow rate can decrease. In the present 
case ,it can he said that density variation is not predominant 
as the flow rate is increasing. Also, to accomodate more 
pressure drop the pressure at the exit decreases. 

(iii) Internal energy was increased hy 15 % at the inlet 

of the secondary circuit ; fig. 5* 3(a) to 5.3(h) will 
represent this situation. 

Ihis increase in energy content of secondary side will 
decrease the heat transfer rate from primary to secondary liquid. 
But due to such a high amount of extra energy on the secondary 
side, quality will increase at the exit. If we go throu^ the 
simple heat balance. 

Heat Transfer Rate (Q.) «= Mass flow rate x Enthalpy change 

% 

With the increase in internal energy the enthalpy at the exit 
is increasing as the enthalpy at inlet is also increasing. But 
there is not much change in H. So the reduction in heat 
transfer rate will reduce the flow rate. In order to reduce the 
pressure drop, pressure at the exit will increase and hence the 
oorr ©sp onding saturation temperature will increase. 
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Coming to tlio primary side, as less heat is transferring 
to the secondary liquid, the internal energy and temperature 
of the outcoming liquid will increase. This increase in 
temperature can slightly reduce the density. Because of the 
isothermal flow assumption on primary side the pressure at the 
exit will decrease correspondingly. This decrease in pressure 
will slightly increase the pressure drop and hence the flow rate 
But all these changes in. primary side are very small compared to 
the secondary aide changes, 

(iv) Mass f low rate was increased hy 20 */. and kept 

constant at the entrance of secondary side (fig. 5.4(a) 
to 5* 4(f)) .will represent this situation. '^ 

As we supply a constant amount of extra cold water this 
reduces the quO'lity on the secondary side. As the quality 
reduces, the flow rate will also reduce. This will result in a 
lesser pressure drop and henoe the pressure at the exit of the 
secondary circuit will increase, increasing the corresponding 
saturation temperature. 

5.1.5 Comparison of Results with Previous Results; 

As there are no experimental results available to 
compare the present results, a trial v/as made to compare these 
results with; thosb of Silva [ 9 ] , But there is" a ma^^^ 
between these two as far as initiai values are concerned. lii 
Silva » a w;ork the transients sra studied f or i the stea^ 
cohditionSj consisting of a total flow rate of about 2^ times the 
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flovr rate we are using. This means, they are considering a 
maximum quality of oray about 0.12 i.e. at the exit of 
secondary of evaporator,. As it is well known, the annvilar 
flow starts at a quality of about 0,2, hence the above region 
is much below film boiling region. So, it oan be seen, that 
always at all the control volumes there will be only subcooled 
(nucleate) boiling. But the major aim of this work is to study 
the transient with two phase flow. So wo have used a quality 
such that, there will be at least 40 ’/. of the secondary of 
evaporator section, in which partial film boiling (annular 
flow) exists. So even if we compare these results with the 
present results we can not see a good agreement. So, even for the 
jpr (33 ent work one case vras studied with the same conditions as 
above for comparison sake and It was found that the results 


were in good agreement. 


This program is not working efficiently for pressure 
variations. If a pressure reduction of 50 % is given at the 
inlet of say primary side then the flow situation can behave 


in two ways. 

At steady state the 
pressures are as shown in. fig. 541 , 1 . 

The pressure drop is around H*/. of 

inlet. Hence jf all of a sudden we 

reduce the pressure to such a low Big.S.l,! 

value, then there can be flow reversal. 



2250 Psi 


270 Pei 


2500 Psi 


But as the present progr{^ is not writ tea to take c^e of the 
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flow reversals it is hard to expect any reasonable results 
from this program for such a hi^ pressure fluctuations. Also 
at the inlet of primary side the liquid is almost at saturation 
point. If suddenly the pressure reduces, there will be two 
phase flov 7 in the primary circuit itself. But as it was 
assumed that always there would be a single phase flow, so 
even a small reduction of about 5 to 10 % in pressure will 
create two phase flow in primary circuit and hence this program 
can not give better results for pressure fluctuations, 

Coming to the secondary side, the reduction in pressure 
at inlet may reverse the flow in the preheater. If all of a 
sudden the pressure at the exit reduces, to adjust to the hi^er 
pressure drop, there can be two phase flow in preheater section. 
But, even here it was assumed that there vrould only be single 
phase flov; in the preheater. So, it is contradicting this 
assumption. Hence, this program can not give accurate results 
for pressure f luctusations on the secondary side. 

Program was run at different step increment on both 
primary and secondary sides. It was found that the progrEjm 
can not be used for hi^^er step increments on primary side. 

On primary side maximum of 10 y. to 12 /. increment in internal 
energy and flow rate can be used. But on the secondary side 
it can give better results oven for very high step increments 
as about 20 % in internal energy and as hi^^i as 30 y ' for 


flow rate. 
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5»1*4 Condusiong g 

Proffl the above discussion and from the graphs it can 
be concluded that the fluctuations in the primary side inlet 
will affect the operating conditions more compared to the 
fluctuations in the secondary inlet flow i.e. the feed water 
flow. Also, the most predominant and important parameter that 
can highly affect the operating conditions, is the primary internal 
energy. As it can be seen from graphs for 5 */. step change, 
maximum rate of variation of exit parameters on secondary side 
is with the 5 */. step change in primary internal energy. Next 
comes the mass flow rate. ' 

On the secondary side, 5 */. step change in flow rate 
does not affect the operating conditions much. This is true, 
because the feed water rate is only 30 */ as compared to 70 */. 
re-circulation. So the net effect in the secondary of evaporator 
is very less. 

It was also found that the mixture parameter's at the 
inlet of the secondary of evaporator are very important in 
finding the operating parameters in a transient case, especially, 
the mixture pressure which will be used as inlet pressure of 
secondary of evaporator'. The interpolated values of derivatives 
also will play very important role in determining the transients. 

Time Step Size s The program was run for 5 different time 
steps 0.01, 0.003, 0.001, 0.0005, 0.0001. It was found that the 
program is hi^ly diverging for the first time step (0,01), 
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A SIMPLIFIED MODEL OF STEAM GENERATORS 
FOR DIGITAL SIMULATION 
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Absiraci — Tliis paper develops a simplified model of a PWR steam generator, A computer programme 
for the steady slate operation was developed, which will be useful for the dynamic analysis for describing 
accident situations. The model incorporates all the various flow regimes and heat transfer regimes that 
arc likely to be encountered by the secondary flow of the steam generator. The primary flow is con- 
sidered as single phase compressed liquid. Given the heat transfer area, pitch and ihe size of the tubes 
the computer programme matches the total power generated within five percent accuracy. Detailed 
pressure and temperature distributions along the length of the preheater and evaporator are also 
computed. 


INTRODUCTION 

With Ihe present interest in the safely studies of 
nuclear reactor behaviour, it becomes important to 
simulate the behaviour of various reactor com- 
ponents. Steam generator forms one of the important 
components of a reactor system where heat is 
exchanged between the primary coolant of the reactor 
and the secondary stream of the working fluid for 
generating power. In the literature there are no 
simplified models available for simulating the steam 
generators excepting the complicated models of the 
reactor thermal hydraulic computer codes such as 
RELAP-4 (Moore and Rettig, 1975) and FLASH 
(Redfleld and Murphy, 1969), etc. Recently few math- 
ematical models were developed by Arwood and 
Kerlin (1977) and Hbkl (1978) to simulate the U-tube 
steam generator. However, these papers use linearized 
equations to predict the behaviour around a steady 
state operating point In reference Arwood and KLerlin 
(1977) the detailed heat transfer regimes are also not 
considered explicitly. The present paper deals with a 
simplified model of the steam generator taking into 
account the various heat transfer regimes and the 
respective correlations into account. A computer pro- 
gramme was developed to analyse a typical steam 
generator utilized in a 600 MW(e) PWR reactor. 

MODEL OF THE STEAM GENERATOR 

The steam generator model used in the present 
study is shown in Fig. 1. The steam generator is repre- 


* On leave of absence from: Indian Institute of Teciv 
nology, Kanpur, India. 


sen ted in two sections, preheater (economizer) and an 
evaporator. In reality the steam generators of the 
PWR reactors have tf-tube geometry, the present 
model simplifies it to parallel tube geometry but 
retains ail (he significant factors such as (he total heat 
transfer area, total amount of heat transferred, pre- 
heater, Lack of detailed data about the internals of 
the steam generator, such as the pressure drop due to 
baffles, quality of steam at the outlet of the generator 
before the steam separators, etc, is circumvented by 
(he pressure balance and heat balance equations as 
described in the later sections. 



Coldstream Feedwater 


Fig. (. Simulation model of steam generator. 
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Table 1 


Mode Flow regime 

1. Single phase 
compressed liquid 

2. Single phase 
compressed liquid 

3. Bubbly flow 

4. Annular flow 


Heat transfer regime 
Pure convection Nu =» 

Local boiling Tp ~ T,, 


Heat transfer correlation 
0.023 Re^-^ Pr^ ** 

0.072 

” exp (p/ 1260) 


Nucleate boiling 

Evaporation by /, = 5(0.00122) A -f- F(0.023)Rc?‘'’ 
forced convection 


A - 




F and S obtained from Chen (1963) 


Reference 


Dittus and Boeltcr 
(1930) 

Tong and Weisman 
(1970) 


Chen (1963) 


5. Annular flow 


Stable film boiling hD^ 


/ I _ aA0*®53 

0-0136 Re nim — 


a 

I -a 




Polomikeui/. (1964) 


As indicated in the figure the coolant from the reac- 
tor outlet enters the primary side of the generator as a 
hot stream in the evaporator section and leaves the 
generator from the preheater section as a cold stream 
to enter the reactor. The heal flow from the primary 
side is conducted through the tube walls of the steam 
generator and enters the secondary stream. The figure 
also indicates that the liquid in the secondary stream 
leaving the evaporator section of the generator along 
with the vapour is recirculated. The feedwater enters 
the preheater section and mixes with the secondary 
stream after the preheater. 

The liquid in the primary side of the stream genera- 
tor is assumed to be in the form of compressed liquid 
transferring heat by pure convection process. On the 
secondary side however the liquid enters the genera- 
tor as near saturated liquid with or without local boil- 
ing. It is envisaged that the secondary liquid will 
encounter bubbly flow and annular flow regimes and 
the heat transfer regimes will be pure convection, 
local boiling. The various flow regimes and the heat 
transfer modes and the correlations used in each 
mode are indicated in Table 1. The logic used to de- 
termine the regime of heat transfer is shown in Fig. 2. 

The computer programme developed calculates the 
equivalent lengths of preheater and evaporator, pres- 
sure and temperature variations along the length in 
both the primary and secondary loops and also the 
exit quality of steam and its variation on the second- 
ary side of the steam generator. The programme 
requires the following input information; geometrical 
characteristics of the steam generator such as the 
pitch, the total heat transfer area, etc,, the thermo- 
dynamic properties of the primary liquid entering the 
steam generator, the total heat received by the 


primary liquid in the reactor (or indirectly the output 
power of the nuclear plant) and the desired vapour 
output and pressure, 

ANALYTICAL MODICL 

In ihc analysis it is assumed that the flow is verti- 
cal. There is no transfer of work acro.ss each control 
volumes indicated in Fig, 1, Further the area of flow 
cross section is assumed to be constant. Each control 
volume is bounded by two junctions, one at the inlet 
and the other at the outlet connecting the subsequent 
control volume. Thus the conservation equation of 
mass, momentum, energy given below representing 
steady state are solved for each control volume: 

Mass: 




Momentum: 


d[pV) dp dF^ ^ fu/v\ (\ 


Energy; 
dt 


A^(pt) = q 


„ 0^ 
dz ' dz 




where q is the heal flux, VF == pAF is the mass flow 
rate, c =« u 4- V^/2 -I- gz is the energy per unit mass of 
the fluid and is the frictional force. A and Ar i‘ep- 
rcsent the flow cross sectional area and the heal 
transfer area, respectively. The above equations are 
cast in the form of volume average quantities and 
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the values at the junction by integrating between the 
functions and defining the volume averages as the 
arithmetic means of the function values. 

The various flow regimes and the respective heal 
transfer correlations are indicated in Table 1. The 
pressure drop calculation required in the momentum 
equation are divided into three categories of friction, 
acceleration and elevation. In reality there are baffles 
placed in the flow of the secondary fluid and thus 
increases the pressure drop due to friction. However, 
since the details were not available, it is proposed to 
use the pressure drop due to friction in the control 
volume //, Apni in the form 


Ap„ = (/ + K) 


2a<, 


Az 

'd/ 



where/,, is the friction factor and X is a form factor to 
account for baffles, change of cross sectional flow 


area, etc., AZ is the height of the control volume and 
Dp is the hydraulic diameter. The above equation is 
valid only for single phase flow. Since the secondary 
fluid will encounter two phases, this equation has to 
be corrected. In this work empirical correlation of 
Baroezy (1966) was used to account for the two-phase 
flow. The work of Baroezy is considered to be con- 
tinuation of Martinelli and Nelson. The empirical 
factor defined as, 

2 ^ (Ap/AL) two phases 
(Ap/AL) single phase * 

was obtained by polynomial curve fitting to the 
graphs around a mass flux of 10® Ib/ft^ h (4.89 x 10® 
kg/m^ h). Thus the friction factor in equation (2) will 
be modified as, 

(/„) two phase = <j>^{f) single phase. (4) 
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The pressure drop due to acceleration was con- 
sidered by the following equation, 

Af.cc = ~(t^n^l - (5) 

dc 

This pressure drop is eAperienced in going from con- 
trol volumes n — 1 to junction ?i. The above equation 
is applicable for the single phase flow. However, in 
two-phase flow slip ratio comes into play. Normally 
the correction factors are dependent on the quality x 
and void fraction a. For annular flow the correction 
factor r is computed as 


L(l-a) a\PvJjill \PvJj 


( 6 ) 


The above equation (5) is modified as 


APacc = — (r„+ iii„ + 1 - (7) 

Oc 


The void fractions are calculated in the case of evap- 
oration by forced convection and annular flow using 
Thom*s (10) correlation, 


l+X(v-l)’ 


(8) 


where y is an empirical constant tabulated as a 
function of pressure. Similarly in the case of stable film 
boiling and annular flow Polomik et aPs (1964) 
correlation is used : 



The elevation pressure drop is simply calculated at 
each junction as 

APdevailon == (10) 

Conduction across the primary tube is calculated 
using the following equation, 

^ r.\n(rjr,) 

where subscripts i and e stand for internal and exter- 
nal while pi and pe stand for primary inlet and exit 
conditions. Kp is the thermal conductivity of the wall 
material, 

METHOD OF COMPUTATION 

Essentially the computer code used in this mode! 
solves the three basic conservation equations of mass, 


momentum and energy in each volume of the primary 
and the secondary side. Conservation of mass can be 
written as — constant. The conservation of energy 
equation can be rewritten as 

I -aS-l''---''-’ 



where is the average amount of heat transferred in 
control volume ii, whose boundaries are described by 
the Junctions n and + 1. Also, H ^ U p/p and 
G = W/A relationships are utilized in the above 
equation. 

Similarly the momentum equation can be simplified 
assuming that in each control volume an average 
value of friction, based on the average flow properties, 
can be calculated. The average liquid properties in 
each control volume are assumed to be the same at 
the interface junction of the subsequent control 
volume. Based on these assumptions the momentum 
equation can be written as 

Pfi ^ Ptf-l "h Pfc.fi) 



K-i)- 


9 

144^;; 


iVn - Vp] 


9 Az 
144^; ' 


The above equation can be rewritten in terms of 
pressure drops due to friction, acceleration and eleva- 
tion as, 

Pfl “ Pff-l “ TC(^P)ft-l (^P)n!]frlcllon 
^PcltvrtllDn* 

The overscores indicate average values and those that 
do not have an overscore indicate the Junction values. 

With the known primary inlet conditions and an 
initial guess at the quality at the outlet of the evapor- 
ator the recirculation rate and an approximate linear 
pressure profile is established, With these one can 
march in the forward direction to determine the 
secondary outlet conditions. If the secondary outlet 
conditions do not match the operating conditions 
furnished, the quality is adjusted and calculations 
repeated until the specified number of iterations are 
exceeded or the total heat transfer and the secondary 
outlet conditions converge within desired limits. 
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ig. 3. Flow chart of the computer programme simulating the steam generator. 


K = 0.7 has given quick convergence for most of the 
results and it is used throughout this work. 

The quality at the exit is taken as 90% of the maxi- 
mum expected value. 

( 12 ) 

When once quality is fixed the secondary flow rate is 
calculated using 

W 

(13) 

^ e,ev 

The recirculation water and feedwater are assumed to 
mix at the entrance of the evaporator at the same 
pressure and temperature. Initial value of pressure 
drop is calculated by making the following assump- 
tions. 

1. Thermodynamic properties are calculated based 
on the exit pressure of the evaporator. 


The major problem in evolving a computer pro- 
gramme is to have the thermodynamic properties of 
water and steam in the various regimes that are likely 
to be encountered in the programme. The present 
programme contains various subroutines in the form 
of polynomials obtained by least square fit to the 
tabulated properties. Similarly some of the two phase 
pressure drop calculations which were available in 
graphical form were also written as subroutines in the 
form of polynomials obtained by curve fitting. The 
flow chart of the computer programme is given in 
Fig. 3. 

The computer programme needs the initial values 
for computation details of which are given below. The 
evaporator height is taken as a fraction of the total 
generator height 

area of heat transfer 

number of tubes x perimeter of each tube ' 

(H) 
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2. An average specific voluine is taken as 


Vs,ev == 




3, The vapour quality is assumed to vary linearly 
along the generator as 


X 


s,ev 



(15) 


Pressure drop in the secondary flow is calculated by 
knowing the flow rates 


P 


3,ev 


1 Gje, 
144 g. 



(16) 


Initially an approximate value is used for pp t>tit 
subsequently it is corrected by using the J value along 
with the corrected thermodynamic properties of water 
at the average pressure in the evaporator* viz, 

hev - (n + iA/))5..v. (17) 

Assuming adiabatic compression, it is possible to cal- 
culate the enthalpies with which the recirculation flow 

and make up water are mixed at the entrance, 

• 


READJUSTMENT OF THE CALCULATED 
VALUi:S (ITERATION) 


After each step the total heat transferred in the 
evaporator is compared with the desired value. If dih 
ference exists beyond tlic desired accuracy then the 
following adjustments arc made. 

(a) Height of the evaporator 

The amount of heat transferred is directly propor- 
tional to the heat transfer area and hence to the 
length of the tube of a fixed diameter. Hence the 
length is corrected as follows; 

fj \ ^ I Q total ^ Qp A 

\UcyJcorrecic«l fr\ \ ‘ 

lota 1 1 calculated ^ ^PA 


(b) Quality of the vapour at the exit of the evaporator 

The total heat transferred is calculated using the 
expression 

e™, - - /„ + - >f) 

“h kPvHHoiir(^^c ^U)s,PA ’ (^4) 

To express Q as a function of quality X, use is made 
of the conservation equation of momentum: 


The amount of heal transferred in the preheater and 
evaporator are calculated using energy conservation 
equations. 


Qm ^ imhe - hdls.PA, (19) 

r 1 1 

Q.V ■= - h,) + ^(vl - (20) 

and 

Qtolal ^ QpA “h Qev* (^1) 

If the total amount of heal transfer differs from the 
desired value, repeat the calculations by changing 
the quality of the vapour at the exit according to the 
equation 

f x,vl + (X, - 2)vl - 2(X„ - l)«,.,t>,,A 

+ {Xc - 4)1),^/- 2(X, - 2)v,,,ii„, Js.e„ 

(2(oial Q>FA kPviipour(^^u,e “ 

The process is repeated till either the desired accuracy 
is achieved or total number of iterations are exceeded. 


IY2 

vaiiQiir 


2A^Xl 


rr(^.- 

{? L 


Vl) 


^Qcji 


^ 2 t/ Apfri 5||(,„)s^ c 


(25) 


It is also assumed that 

(^ti>)s, preheater ~ (^t /,(»)?, evaporator* 
iVe) a,|Uc|tcaCcr (i»i..) s.evaporntor* 

Then the change in total heat flux dQ is set equal to 
zero to obtain 

2^0X11 

f Apfriciot)/ , ^ 

1 — 

Ap , 

A u V A u / I S, ey 




(26) 


(c) Pressure, temperature and enthalpy of the liquid 
entering the evaporator 

Pressure drop in the evaporator is the sum of the 
pressure drops in each control volume (calculated by 
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Tabic 2, Characteristics of the U-lnbe steam generator 


Number of lubes 

External diameter of the tube 

Nominal thickness of the tube 

Total area of heat transfer 

Tube spacing 

Pitch 

Heal transferred 

PRIMARY 
Pressure at entrance 
Inlet temperature 
Coolant flow rale 

SECONDARY 
Make up water pressure 
Make up water temperature 
Vapour pressure 
Vapour flow rate 


4674 

0.75 in, {1.905 cm) 

0.043 in. (0.109 cm) 
48,300 ft^ (4487 m"') 
1.0675 in. (2.711cm) 
Square 

3.211 X lO^Btuh"^ 
(0.941 X lO'^W) 

2250 psi (155.27 kgf cm^) 
620"F (326.7X) 

3.55 X lOMbmh”' 

(1.61 X lO^gh”^) 

1121 psi (77.36 kgfcm-2) 
430"F (221. rC) 

920 psi (63.49 kgf cm 
4.08 X 10^1bmh“* 

(1.85 X lOngh"^) 


v*''Or'joor*^oo-mO>^ONro 

c>doodoc>odo 


»-iri»'-roooociOoo^O'oor^ 
^ ^ Tf Tt* ^ 


.^»-^xooooooNor^o>oo ^ 
^'drv:odoS\dr^‘ri»o'oK 


OMt^r'ir^ONOOfnr-inin 
,, M m ro fo Tf en m fO rr> ri 

X fvj (vj tsj (vj r-j Cl 


the programme). Using this pressure drop the 
enthalpy of the liquid entering the evaporator can be 
calculated. Knowing the pressure and enthalpy the 
temperature is calculated. 


gooooit-ri^Tfoon^-'O 

pS; 

> aH oi ri r-i <n ni H <n ri 


(d) Heat transferred in the preheater 

Knowing the new value of the enthalpy of the 
liquid at the entrance of the evaporator the amount of 
heat transferred in the preheater is recalculated using 
equation (19). 

(e) Secondary flow in the evaporator 

Knowing the corrected quality, the flow rate of the 
secondary steam In the evaporator can be calculated 
using equation (13). 

The above corrections are repeated till the desired 
accuracy is achieved in the calculation of the total 
heat transferred. 


The above programme is utilized for analysing the 
steam generator of a typical PWR plant whose data is 
given in Table 2. 

Results of the computation are given in Table 3 and 
the resulting errors of the various parameters are 
listed in Table 4. From the table it is clear that the 
programme calculates most of the design parameters 
with less than 10% error, The error in the total heat 
flux is less than 1% and the length and hence the heat 
transfer area are reproduced within 6% error. Also it 
should be pointed out that in most of the cases con- 
vergence was achieved in less than eight iterations, A 
typical run took less than 1 min on an IBM 370/145 


|2+, 

E X *0 




OOCpOOWOOOOt30 009 C<» 

Sj Jo cJoodooood 




cjs 'O — < 4^ r^, 00 f<> 'St 

T d d H d d ^ 


fni^oothor^f^'oOoJ. 

oqo&io 

‘‘Jvdvodvosdd'Od'Ovo 


^sqvdvor^r^r^oooocooq 


.0*/^00'^0ow^0«r» 
rN «N cn rs m 
dddddddddd 


«<Ncn^»Ovor^OOOvO 
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Table 4. Comparison of the results of steam generator simulation 


Maximum percentage 

Value deviation 


Parameter 

Expected 

Average 

Minimum 

Maximum 

Expected 

value 

Average 

value 

i-rA 



16.200 

16.174 

16,241 

— 

0.25 


— 

34.960 

33.709 

37.009 

— 

5.86 

Apr 

-- 

18.25 

17.70 

19,10 


4.66 

^Ps, ee 

— 

3.95 

3.87 

4.30 

— 

8.86 

Qpa X 10-® 


0.49843 

0,49843 

0.49843 

— 

0.002 

Qc. X 10-’ 

— 

2.7401 

2.7308 

2.7445 

— 

0.34 

Qm X 10 

3.2115 

3.2318 

3.2992 

3.2429 

0.98 

0.34 


48.300 

46.951 

45.79S 

48.821 

5.18 

3.98 

Wy X 10'* 

4.085 

4.1141 

4.1062 

4.1268 

1,02 

0.31 

X. 

— 

0.2908 

0.1979 

0.3834 


32.00 


computer. Several test runs were made by varying the 
pitch and spacing and the above trend persisted in all 
the runs. 

The interaction of the primary and secondary cir- 
cuits will be crucial in transients when the steam 
generator will be operating outside the range of 
steady state design. This programme is being 
extended in order to study the transients in the 
primary circuit, such as reduced mass flow rates and 
temperature perturbations of a PWR type nuclear 
power plant. 
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APPENDIX II 

This consists of description of some important sub- 
roMtines used for transient simulation. 

(1) CP AGO s It calculates specific heat Cp of ccmpressed 
water at constant pressure as a function of temperature 
(°P) at various pressure ranges. The unit of specific 
heat is Btu/lb~°F. 

Cp = A.T^ + B.T' + C 

Values of A,B,C for various pressure ranges are tabulated 
in Table 1. 

(2) CP AGS 5 It calculates specific heat Cp of saturated 
water as a function of temperature T (*^P) 

Cp = A.T^ + B.T + C 

The values of constants A,B,0 for various temperature 
ranges are tabulated in Table 2. 

(5) CTAGC s It calculates the conductivity (K) of compressed 
water in Btu/hr~ft~°E as a function of pressure (P si) and 
temperature (^) 

K = A.T^ + B.T + C 

The values of constants A, B,C for various pressure ranges 
are tabulated in Table 5* 
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(4) CTAGS ; It calculates the conductivity (K) of saturated 
water as a function of temperature T (^) 

K = A.T^ + B.T + G 

The values of constants A,B,C for various temperature 
ranges are tabulated in Table 4. 

(5) G5MTP 2 It calculates the conductivity of the material 
of the tube. Theinaterial chosen in this problem is 
Inconel-718. 

K = 57.78 (-1.98x10“^ + 8 . 91x10“'^ + 0.151) 

where K is in Btu/hr-ft-^P and T in °C. 

(6) PTMDP ; It calculates the multiplication factor needed 
to evaluate t\^ro->phase frictional fpressure drop. P or this 
Baroczy's plot [11] is used. The x axis of this log- log 
plot is divided into 6 regions and the segments in each 
region is assumed linear such that 

where A»B are constants 

Region I ; 0,0001 I 0.001 

Region II s 0.001 1 0,003 



7.7 


Region III 

0.0u3 

I 

0.007 

Region IV 2 

0.007 

I 

0.05 

Region V j 

0.05 

I 

0.1 

Region VI 

0.1 

I 

1.0 


(7) QJFOCV ; In calculatea P , S, and cr to be used in 
Chen's correlation 

I = Exp [a + b.ln (x;^^)3 

Values oi a,b for different ranges of are tabulated 

in Table 3 (a) » 

S = c . + d . B + e 

vfhere, B = Re^ . j^l*25 

Values of Cj d, e for different ra;:iges of B are tabulated 
in Table 5(b), 

2 ■ ■ 

cf = a . P + b.P + c 

Values of constants a, b, c are tabulated in Table 5(c). 

(8) TLSUB i It determines the temperature of compressed 
water for a given pressure (P, P si) and specific enthalpy 
(H, Btu/lbra ) . 

H = A.P^ + B.P + 0 

Values of the constants A, B,C for different temperatures 
are tabulated in Table 6. 
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(9) VIAGC ; It determines the viscosity ([a) of compressed 
water in. Ibm/hr-ft as a function of pressure (P, psi) 
and temperature (T, *^) 

li = A + B + C T + D 

Values of constants A» B, C, and D for different pressure 
and temperature ranges are tabulated in Table 7. 

(10) VIAGS ; It determines the viscosity of saturated wafer 

^ as a function of tanperature (T) 

11 = A + B T^ + 0 T + D 

Values of constants A, B, C ar^d D for different temperature 
ranges are tabulated in Table 8. 

(11) VIVAS i It determines the viscosity of saturated^^^^^^^^^^^^^^ 
vapour in Ibm/hr-ft as a function of pressure (psi) 

H = 1.15826(-0. 24553.10“® P^ + o;i3778. 10“^ P + 0.02910^9^^ 

(12) DERIV 5 This subroutine evaluates the various derivatives 
of pressure and temperature with respect to internal energy 
and density at constant density and internal energy 
respectively. These interpolated values are tabulated in 
two tables which are Table 8( a) and Table 8(b) . 
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APPKNDIX UT 


CB^:CK1N(5 'Vm PPOGHAM 

IMPUT SPKCXFiCATICiWS OF THE STEAM (iP^NERATOK 
-.^flJMIJER OP’ COWTFUU- VOLUMES 10 

"PRECISION HEOUIRPriD O.IOOOOE-Ob 

-MAXIMUM WO, OF ITEIRATIOWS 25 

-ACCELARATIQN DUP: ID GRAVITY U,31174EFU2 

PRINT OPTION 0 


CHARACTPIRISTICS OF THE STEAM GENERATOR 


NUMBER OF TUBPJS 4674 

EXT, DIAMETER 0,75000E+00 

INT, DIAMETPIR 0.66400p:+00 

TOTAL HEAT TRANSFER ARP:A 0,48300E+05 

FORM COEfF, 0,27144E+01 


SPACING 0,10675E:+01 
RECTANGULAR TUBES 


PRIMARY CIRCUIT 


PRESSURE AT entry 0,22500E+04 

TEMP, AT E.NTRY 0,6l830E+03 

ENTHALPY AT ENTRY 0,64040E:+03 

FLOW RATE 0,35b50E+08 

fi::STIMATKD STEAM PRESSURE 0,22300E+04 
SECONDARY CIRCUIT 


PRPJSSURE 0,11210P>04 
TEMPERATURE ().42a95E+03 
PJWTHALPy 0.40780E+03 
VAPOUR PRESSURE 0.92000Et03 


gUALlTY 0,31900E+00 

VAPOUR FLOW RATp: 0,40a50E+07 

KATE OF HEAT TRAN5Ffe:R 0,32H5K+10 





TABLE 8(a) 


♦ SPECIFIC INTERNAL ♦ OERIi/ATIVE OF * DERXVATXVE OF » 

♦ ENERGT » PRESSURE * PRESSURE ♦ 

i (U) * (dP/dP) * (dT/dP) * 
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♦ 

♦ 

♦ 

♦ 

* 

* 

♦ 

♦ 

♦ 

>»< 

♦ 

* 


68,04 

168,04 

269,52 

347,27 

428.60 

485.10 
508,50 

532.60 
557,40 

583.10 
609,90 

638.30 

668.70 

702.30 

741.70 


♦ 

♦ 

» 

♦ 

* 

♦ 

♦ 

* 

* 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 


5577,95 

5748,97 

5186,20 

4230.15 

3656,88 

2965.79 

2712,02 

2466,91 

2178.00 
1949,84 

1667.01 
1428,93 
1201,41 

972,17 

999,50 


♦ 

♦ 

* 

4 

♦ 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 


1.89075 

6*^0073 
8,61391 
9.37016 
10.02420 
10,21480 
10,50280 
10,56390 
10,81580 
10,71360 
10,88470 
11,02110 
11,08470 
13, 66450 


4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

* 

4 

4 


^4]\f:^;t4444 444444444444444444444444444444444‘4444444444444444*4444444 


TABLE 8(b) 


tlp/du. s A + B , U 
dT/du a C + 0 . U 


DENSITY 
( ^ ) 

61,9963 
60,1178 
57.3132 
53,6538 
51,4569 
48,9476 
47,8240 
46,6200 
45,3104 
42,3191 
40,5681 
38.5654 
36,1402 
32,9815 


444444444444444444444444444444444444444444*44444444444444 

4 A * B ♦ C 4 D 

4 4 4 ■ ■ ■ »'■" ^ .'^■■ 4 . 

:^^44444444444444444444444444444444 4444444444 44444444444444 



¥ 

81.1512 

¥ 

¥ 

0,125973 

¥ 

¥ 

0,978686 

¥ 

7 

,99847E-4 * 

* 

175,7350 

¥ 

-0,167236 

¥ 

0,941878 

¥ 

7 

.99847E-4 4 


296,9620 

¥ 

-0,460533 

¥ 

0.987761 


7 

,67052E-4 * 


517,0680 

¥ 

-0,915354 

¥ 

1,271330 

¥ 

2 

,99893E-r5 ♦ 


593.3069 

¥ 

0,993442 

¥ 

0,635391 


4 

.44417E-3 ’If 


3.54932 

¥ 

0 ,284503 

¥ 

0,286995 

¥ 

2 

,1309tE-3 f 

¥ 

72,99120 

¥ 

0,112449 

¥ 

1,148570 

¥ 

3 

>86180E-4 ♦ 

¥ 

-65,34110 

¥ 

0.352303 

¥ 

1,374650 

¥ 

-2 

,57190E-S f 

¥ 

-152,7980 

-339.2240 

¥ 

0.477300 

¥ 

0,572154 

¥ 

1 

,40439E-3 * 

¥ 

¥ 

0,710325 

¥ 

-2,20264 

¥ 

5 

,77373E-3 * 

¥ 

-70,2542 

¥ 

0.249500 

¥ 

0,383808 

^ ;¥:■ 

i 

,5l964E-3 » 

¥ 

767,6040 

¥ 

-1,014860 

¥ 

4,878320 

¥ 

-5 

.19858E-3 ♦ 

¥ 

74,7736 

¥ 

0.000000 

¥ 

22,23570 

¥ 

-0 

,5294845 * 

¥ 

79,1472 

¥ 

0,000000 

¥ 

1.67739 

. ¥ 

0 

.ooooooor 

¥ 


¥ 


¥ 


¥ 


-■ir ' ''■. 4 ‘ 



NtfMBEFl 0&' Xl’KHATiaNS 21 
PKECISIOi^ OftTAlNED 0.75287E-O7 

£QlJIVAt.eHJT CiENGTH UF PREHKATKR 0,15984E+02 
TOTAIi EQUIVALENT LENGTH OF' EVAPORATOR 0,34786E+02 
TOTAL EQUIVALENT LENGTH 0,50770^+02 
TOTAL HEAT TBASFER AREA IN SQ.FT. 0,4b593E+05 

DlwSTRiaUTION OF PRESURRE ANU TEMPERATURE IN THE PRIMARY CIRCUIT 


PRESSURE FLU ID TEMP 

CPSI) COEGREE F) 

0.22500E+04 0,61R30E+03 

VOLUME .1 a,22492E + 04 0,6l389E+03 

0,22485E+04 0,60948E+03 

VOLUftF. 2 0,22477B + 04 0.60557E+03 

0.22470E+04 0,60165E+03 

VOLUME 3 0,22462E+04 0.59776E+03 

0.22455E+04 O.59387E+03 

VOLUME 4 0,22447E+04 0.59035E+03 

0,224,40E+04 0,5B684E+03 

VOLUME 5 0,22432E+04 0,58379E+Q3 

0,22425E+04 0.58073E+03 

VOLUME 6 0,2241RFJ+04 0.57808E+03 

0,22411(2+04 0,57542E+03 

VOLUME 7 0,22403E+04 0,57321E+03 

0,22396E+04 0.57100Et03 

VOl.UME 8 0,22389K+04 0,b69l0E+03 

0.223U2E+04 0,56720E+03 

volume 9 0.22375E+04 0.b655iE+03 

U.22368E+04 0,5b387E+03 

VaLUME 10 0.22360E+04 0.56239E+03 

0,22353E+04 0,56092E+03 

PREHEATER f),22337E + 04 0,55535B:+03 

0,2232.1E+04 0,S4979E+03 


WALL TEMP, ENTHALPY FLUX 

COEGREE F) (8TU/LB) (BTU/SQ.FT/HR) 

0.64040E+03 


0,58949E+03 

0.63383E+03 

0.62727E+03 

0.14625E+06 

0,58389E+03 

0.62146E+03 

0,61566E+03. 

0.12923E+06 

0.5785eB+03 

0.61058E+03 

0,60549B+03 

0,1 T324E+06 

0.57358E+03 

0,60i07E+03 

0,59665E+03 

Or9a4i2E+05 

0,569i5E+03 

0.59282E+03 

0.588g8E+03 

0r85479E+05 

0.56529E+03 

0.58564B+03 

0,S8230E+03 

O,7 015E+O5 

0.56247E+03 

0,67951E+03 

0.57673E+03 

0.61;947E+05 

0,55988E+03 

0,57435E+03 

0.57196E+03 

0.53182e+05 

0,55746E+03 

0,b6987E+03 

0,465786+05 


0,56777E+03 


0.55522R+03 

0,56592S+03 

0,56406E+03 

0,41:359E+ij5 

0,54952E+03 

0,55706E+03 

0.55006E+03 

0,33924E+05 



(z ^ 0 ) - rv) \G, 



